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In the context of cosmological perturbation theory, we derive the second order Boltzmann equation 

describing the evolution of the distribution function of radiation without a specific gauge choice. The 

essential steps in deriving the Boltzmann equation are revisited and extended given this more general 

framework: i) the polarisation of light is incorporated in this formalism by using a tensor-valued 

distribution function; ii) the importance of a choice of the tetrad field to define the local inertial 

frame in the description of the distribution function is emphasized; iii) we perform a separation 

^^ between temperature and spectral distortion, both for the intensity and for polarisation for the 

, ^ first time; iv) the gauge dependence of all perturbed quantities that enter the Boltzmann equation 

^^ is derived, and this enables us to check the correctness of the perturbed Boltzmann equation by 

explicitly showing its gauge-invariance for both intensity and polarization. We finally discuss several 

C implications of the gauge dependence for the observed temperature. 

'Z^ PACS numbers: 98.80 

u 

' The non-Gaussianity in the Cosmic Microwave Background (CMB) has been one of the hottest topics in cosmology 

^ because it could open a new window for probing the primordial universe. Recent CMB observations, especially 

"^ WMAP [T] and Planck P], have confirmed to a very high accuracy that the primordial curvature perturbations have 

Ci a nearly scale invariant initial power spectrum and the associated statistics is nearly Gaussian. These observations 

'"^ are consistent with the predictions of an early inflationary era driven by a single slow-rolling scalar field. 

-v^ The possibility of non-Gaussianity in the primordial curvature perturbations was discussed for the first time quan- 

^>. titatively by Komatsu and Spergel f3]. They parameterized the level of non-Gaussianity in the potential $(a;) (the 

i^\ curvature potential in the Newton or Poisson gauge) by 

(N 

g ^x) = <i>u^) + f^{:^'[<Plix)-{<i>l{x))\ . (1.1) 

'^ Here, ^h{x) denotes the Gaussian part of the perturbation, or in perturbation theory its linear part, and (•■■) 

~r designates the statistical average. This type of non-Gaussianity leads to a non-vanishing three point correlation 

I function, or equivalently in reciprocal space to a non-vanishing bispectrum. The prediction for the possible values of 

• • this parameter /^l^' from a phase of single-field slow-roll inflation was first performed by Maldacena [4] , and it was 
• ^ shown that it is of order of the slow-roll parameters and thus highly suppressed. Hence, if a /^l''' of order unity or 
K> greater is detected, this simplest model of inflation will be ruled out. 

;_( Unfortunately, the interpretation of the measured non-Gaussianity is not so straightforward because we do not 

^ observe directly the primordial non-Gaussianity in the curvature perturbation but its effect on the CMB fluctuations. 
Therefore, to relate the primordial curvature perturbation to CMB, we first have to compute the evolution of pertur- 
bations after inflation. These effects can be split unambiguously in two parts: i) a linear transfer that cannot create 
a non-Gaussian signal if the initial conditions are purely Gaussian, and ii) a non-linear transfer that generates a 
non-Gaussian signal in the observables even if the initial conditions are purely Gaussian. The resulting non-Gaussian 
signal from i) is often called primordial non-Gaussianity and all the possible sources of non-linear evolutions which 
enter the category ii) are called secondary non-Gaussianity. Recently the Planck collaboration provided a constraint 
/nl*^' = 2.7 ± 5.8 [2]. This result was obtained by subtracting one of the secondary non-Gaussianities that arises from 
the correlation between the lensing and integrated Sachs- Wolfe effect. This clearly demonstrates the importance of 
subtracting all the secondary non-Gaussianity consistently in order to obtain an accurate constraint on the primordial 
non-Gaussianity. 

The evolution of the perturbations on super-horizon scales is well understood, even fully non-linearly using either 
a covariant approach [SHE], or a separate universe approach with the so-called 6N formalism [SHU], since it leads to 



a conservation law for the curvature perturbations in the case of adiabatic perturbations. On the other hand, the 
evolution for modes below the horizon scale is not so simple analytically, especially at the non-linear order, and the 
use of a kinetic description cannot be avoided on small scales since radiation starts to develop an anisotropic stress. 
In order to obtain numerical results for the non-linear evolution, we need to derive and solve without approximations 
the coupled system of non-linear a) Einstein equation for the metric, b) conservation and Euler equations for fluids 
and c) Boltzmann equation for radiation (photons and neutrinos). Note that the conservation and Euler equations 
can always be deduced from the lowest moments of the Boltzmann equation, and the full set of equations is often 
only referred to as Einstein- Boltzmann system of equations. 

In order to follow this roadmap, the second order Boltzmann equation was written down in the Poisson gauge in 
Refs. [12141 7j . The gauge dependence of the distribution function was obtained at linear order [18] and then at second 
order [TH] but leaving aside the problem of polarisation. It was then extended to include polarised light in Ref. |14| . 
The system of equations was then solved numerically in Fourier space in Poisson gauge in |16| , and it was reported that 
the secondary effects around the last-scattering surface could mimic a primordial signal of /^l'^' ^ 4. Recently there 
have been a huge progress in improving the numerical calculations and clarifying the amplitude of various secondary 
non-Gaussianities at recombination |20H22| . and a consensus emerged that when including all the non-linear effects 
around recombination and the integrated early effects after recombination, it could mimic a primordial signal of 
/nl**' ~ 0.8, as expected from analytic approximations [23] . 

The description of the spectral dependence of the distribution function is also crucial at second order. Indeed, at first 
order there are no spectral distortions and the perturbation of the photon distribution function can be understood as a 
single, spectrum-independent temperature fluctuation. However, at second order, there appears a deviation from the 
Planck distribution, resulting in a continuum of spectral distortions, which in principle must superimpose the thermal 
Sunyaev-Zeldovich effect [24] . To describe this distortion, we use a direction and position dependent Compton y 
parameter [25[ I26j , and also introduce a similar tensor- valued variable to describe the distortion in the polarisation, 
thus extending the formalism introduced in Ref. [55] . In this paper, we derive the second order Boltzmann equation 
without restricting to a specific gauge, and including polarisation. Reflecting on the above, our motivation is two-fold. 

First, since the structure of the second order Einstein and Boltzmann equations depends very much on a choice of 
the gauge, we have to find a gauge in which we can numerically solve this system accurately and quickly. Therefore, 
it is preferable not to specify the gauge from the beginning but to formulate the equations without specifying it. 
We can impose different gauge restrictions in their final form to explore the stability and efficiency of the numerical 
integration. Second, we would like to check the equations derived in Refs. [1311111113. As a direct check, we recover 
them in the specific case of the Poisson gauge. Then, as an indirect check, we revisit the transformation properties 
of the distribution function and the metric perturbations and confirm that the perturbed Boltzmann equation is 
gauge-invariant up to second order in perturbations, thus increasing our confidence in the rather lengthy derivation. 

In Ref. [22] , it was found that the inclusion or omission of certain line of sight terms can make a large impact on 
the estimation of the bias to the primordial non-Gaussianity due to the secondary non-Gaussianity. In Refs 1201 [55] 
all physical eflFects were included except for lensing and time-delay. These time-integrated effects require a separate 
analysis because at later times small-scale multipoles get excited and numerically it is very difficult to evolve the 
equations. In this paper, we point out that the separation of these effects depends on a gauge. Given that the 
lensing-ISW cross correlation gives the largest bias to the primordial local type non-Gaussianity, one should bear this 
gauge dependence in mind when separating these time integrated effects in the calculations. 

The choice of the gauge and the associated choice of the tetrad field for the distribution function is also crucial in 
the interpretation of the quantities as observables. These subtle details do not affect our interpretation of observables 
in the linear theory since it is only relevant for the monopole and the dipole. However it is no longer the case at 
second order in perturbations. Wc must understand the transformation properties of the distribution function under 
a gauge transformation or a change of the inertial frame and determine what is a gauge and a choice of inertial frame 
that is related to CMB experiments. 

The structure of this paper is as follows. In section [llj we give the definitions of the variables that we use for 
the metric, momentum and distribution function. Especially, to express the perturbation of the metric, we use 
a geometrical (3 -I- 1) decomposition, or the ADM ?I7] parametrisation of the metric. At first order, there is no 
particular advantage in using this formalism, but various expressions are simplified at second order for the choice of 



the inertial frame that we make. In section III we derive the second order Boltzmann equation with polarisation 
without restricting to a specific gauge. In section jTvl we discuss the gauge dependence of the variables. We carefully 
investigate the gauge transformation of the metric, momentum and the distribution function. We then check explicitly 
the gauge invariance of the perturbed Boltzmann equation up to second order as a consistency test. Finally, in 
Section Ivj we summarize our results and we comment briefly on the relevance of our formalism for the observed GMB 
anisotropy. Useful technical details are gathered in the appendices. 



II. DEFINITIONS 

In this section we build all the tools which are used for the description of polarized radiation in cosmology. We 
first review briefly the parametrization of cosmological perturbations, and explain how a photon momentum can be 
uniquely described by its energy and direction once a suitable tetrad choice has been made. We then introduce the 
tensor-valued distribution function which is used to treat statistically a gas of polarized photons, and which is the key 
object in the Boltzmann equation, and we finally present how it can be decomposed into its main spectral components. 



A. Spacetime coordinates and local inertial frame 



We shall use the ADM formalism to write down the expression of the perturbed metric where the metric can be 
decomposed as 






(2.1) 



where N is the lapse function, /3* is the shift vector, ^ij is the spatial metric, and indices of the Latin type (i, j, k ■ ■ ■) 
run from 1 to 3. To describe the perturbations around the flat Friedmann-Lemaitre-Robertson-Walker (FLRW) 
space-time, perturbation variables, a and /i^j, are introduced as 



N =1 + a, lij = 6ij + 2h.ij . 
For simplicity, we use the following definitions 



(2.2) 



(2.3) 



where the spatial indices are raised and lowered with 5ij and (5'-', rather than with 7,^ and 7*-'. As is clearly seen 
below, the ADM form of the metric perturbation will simplify the expressions of the perturbed Boltzmann equation. 
Any perturbation X will be expanded into its first and second order parts as 



- y(i) 



X = X 
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-X 



(2) 



(2.4) 



The relations between the ADM variables and the usual definitions of cosmological perturbations are provided in 
Appendix {M 

The Boltzmann equation is better formulated by explicitly using a local inertial frame at every point of the space- 
time and this can be achieved by using a tetrad field. It is a set of four vector fields which satisfy 



^(a)(6) = ffM'^eCa) e(b) 



dt^u - V{a){b) 



.{a) Jb) 



(2.5) 



These conditions determine the choice of tetrad only up to rotations and boosts. Here the following particular tetrads 
are chosen up to second order accuracy 



s(o)^ = a{-N, 0, 0, 0) , e^'V ^ a ( f3' + h'jP^,S'j + h 



h^'^hk, 



and the inverse tetrads are given by 



e(o)^ = -e(°)'' = - 



1 / 1 



/3' 



a VVV N 



Hi) 



0,%)(fc)(<5^"-/j-'''+2/i''/i/' 



(2.6) 



(2.7) 



The time-like tetrad is chosen to be orthogonal to the constant time hypersurfaces since e'^^) ex drj. As for the 
spatial tetrads, this choice corresponds to asking that there is no rotation between the background and the perturbed 
tetrads [H]. 



B. Momentum 

To facilitate the separation between the magnitude of the momentum and its direction in a covariant manner, let 
us consider the projection of the momentum of photon p^ onto the set of tetrads, 

We introduce the conformal momentum of photon rather than the physical momentum p^'^' 

g(°) EE apt'') . (2.9) 

Since the momentum of photon satisfies the null condition p^p^ = 0, or equivalently q^°''qia) = Oi only three compo- 
nents among four are independent, that is 

g('')g(,)=0, ^ (9(°))' = %)0)9<^'9<^^ (2.10) 

Thus the three spatial components q^^' can be regarded as such independent variables. Furthermore, q^^^i can be 
decomposed into its magnitude q and direction n^*' as 



^y^%)(,)gWgO-) = 1^(0)1, „W^^. (2.11) 



Physically the above q can be understood as the conformal (re-scaled) energy, q — a£'phys, seen by an observer 
orthogonal t o tim e constant hypersurfaces. 



From Eq. (2.9 1, the components of momentum p^ are expressed as functions of (q, rv-"^') up to the second order as 

r(l-a + a2), (2.12a) 



a" 



( ?i(*) ~P' ^ h'ju^^'^ + a/3* + -/i^'^/ifej-n^^'M . 



p' = ■% ( ?i^*^ - /3' ~ h'ju'^^^ + a/3* + -h'^hkjn^^^ ] . (2.12b) 



Conversely, {q,n^^') are given by the components of momentum as 

q = a^{l + a)p° , (2.13a) 



n"- 



(1 - a -I- a^)5*j + (1 - a)h'j h'^'^kkj 



^ + {l-a)P'+P^h',. (2.13b) 

pU 



One can introduce a projection operator in terms of e^'^'^ and n^. The projection operator, often called the screen 
projector, is defined as 

V = gM- + e(°Ve^°'--"M«-. (2.14) 

where n^ , the direction vector of photon, is defined by 

n^' = e(,)'^n(*) . (2.15) 

Clearly 5^1^ is a projection of the tangent space onto a two dimensional plane orthogonal to both e^^'fj, and n^ 
since S^'^e^'^'^ ^ and S^^n^ vanish. Its expression in tetrad components reduces necessarily to the identity of the 
two-dimensional subspace which is left invariant by the projector, that is 

%)(j) = hi){3) - "-(»)"-0) ' ^{o){o) = '5'(o)(.o = . (2.16) 

C. Distribution function for photons and Stokes parameters 

In order to describe the polarisation of radiation, we introduce a tensor- valued distribution function /^i,, which 
is complex valued and Hermitian. The construction of this distribution function is discussed in Appendix |B] It is 
independent of the choice of the electromagnetic gauge and contains only four physical degrees of freedom since it 
satisfies the conditions 

^.6(0)^ = /^.e(o)'' = /^.n^ = I^^.n- = . (2.17) 



Note that the distribution function depends on the observer's velocity, u'^ = e(o)'^j used in its definition. As long as no 
confusion arises from such dependence, we omit to specify it. In the case where this is needed, mainly when studying 
the transformation properties of such a quantity, we shall use the notation /^i"' to stress that the tensor-valued 
distribution function is dependent on the observer's velocity and thus on the choice of the tetrad field. 

The four degrees of freedom can be extracted by decomposing /^^ into a trace part, a symmetric traceless part and 
an antisymmetric part as 

/^, = hs^, + P^, + '-ep^.ae^ofn'^V , (2.18) 

where the antisymmetric tensor is defined by 

£afijS = <^[afSiS] , £0123 = V^ , Or e(o)(i)(2)(3) = -e(°)(i)(2)(3) ^ ^ ^ ^2.19) 

/ is the intensity and V is the degree of circular polarisation. P^^, encodes the two degrees of linear polarisation 
(so called Q and U Stokes parameters). All these functions, together with the original tensor- valued distribution 
function, are functions of the position on space-time x^ = (77, cc') and on the point in tangent space. This point in the 
tangent space can be chosen to be parametrized either by the components p^^ in the basis canonically associated with 
the coordinates system, or alternatively by their Cartesian counterparts p^'^' . In fact we will choose to parametrize 
the tangent space by the components of the conformal momentum in tetrad space, q^"^' = ap^'^\ expressed in their 
spherical coordinates q and rS^\ as this leads to the most simple form for the Boltzmann equation as we shall see 
further. 

D. Spectral distortion 

On the background space-time, the distribution function, which is characterized only by the intensity /, is given by 
a Planck distribution whose temperature T depends only on 77 due to the symmetries of the FLRW universe. As we 
will check later, the background temperature scales as ex 1/a. We thus have 



-^('7, 9) = -^BB 



a{ri)T{T]) 



with /bb(:c)= ,, ,, . (2.20) 

(e^ — 1) 



At first order in perturbation, the fluctuation of intensity can be described as a fluctuation of temperature ST which 
is independent of q. There are two reasons for this. First, as we shall discuss further, gravitational interactions do 
not induce spectral distortions in the sense that they shift all wavelengths by the same ratio. Second, the collisions 
at linear order in perturbation do not induce spectral distortions and the redistribution of the photon directions 
resulting from it can be described by a direction dependent temperature. A similar procedure can be followed for the 
description of polarisation at first order. 

However at second order the situation becomes more complicated since the Compton scattering at this order of 
perturbation induces spectral distortions which cannot be reabsorbed in a simple direction dependent temperature. 
As a result, the photon distribution is not described by a Planck distribution function, but fortunately it is sufficient 
to use two direction dependent quantities. The first remains the temperature and the second describes the type of 
spectral distortion generated at second order. Actually in general, at the n-th order, n directional dependent functions 
would be needed i2_5i l26j to characterize fully the spectrum. 

In order to parametrize this distortion, we introduce on top of the temperature T, the so-called Compton y pa- 
rameter. In this section we will omit the dependence of all quantities on the coordinates x^ and we will focus on the 
dependence on the tangent space coordinates (q,n(*-'). The distribution function can be expanded around a Planck 
distribution in the so-called Fokker-Planck expansion as [IS] 



d 






^/bb(^)+.(.«)P,^/bb(|.), (2.21) 



where 



-.2 _ „-3 d ( ^^ d \ 92 ^ Q 



'' d\nq \ d\nqj d\nq'^ dlnq 



Because the number density of photon is given by n ex a~^ J Iq^dq, the y term does not contribute to the photon 
number density and the temperature T is the temperature of the black-body that would have the same number density 
(see Ref. [T5] for a discussion on other possible definitions for the temperature) and we call it here number density 
temperature. It can be expanded around the background temperature as 

T(n«) = t(77) [l + e(n(''))l . (2.23) 



Note that the expansion (2.21 ) is not the same as Eq. (11) nor Eq. (15) of Ref. [25 . Indeed, the temperature of the 
Planck spectrum around which we expand is neither the physically motivated logarithmic averaged temperature of 
Ref. [25j nor a fiducial temperature, but another physically motivated temperature (the number density temperature) 
that suits better to describe the spectral distortion of the type that appears in CMB. 

However, when performing perturbations in cosmology, we need to refer to the background space-time temperature 
T, not to the local number density temperature. Thus it is convenient to expand the distribution function around a 



Planck distribution at T rather than T. Expanding Eq. (2.21 ) in Q up to the second order, we obtain the expansion 
as 

/^/=e(i)-{e + e^)^/e=(X) + (.+ ie»)l,^/„e(X). (2.24) 

where we used the fact that y is at least a second order quantity. Here, in order to simplify the notation, it is implied 
that Q and y depend on a;^ and n^'^ For a given /, the spectral components Q and y can be extracted by performing 
different types of integrals on q (see appendix Ie] for details). This expansion is similar to Eq. (11) of Ref. ^25j when 
only second derivatives of the Planck distribution are kept. 

Now we want to obtain a similar decomposition for polarisation. Indeed, when dealing with polarisation we also 



need to expand its spectral dependence in a way similar to what has been performed for the intensity in Eqs. (2.21 ) and 



V 



(|2.24 ) , that is we want to separate the polarisation tensor into a spectral distortion Y^^, and non-distorted component 
However, this separation is slightly different given that there is no polarisation on the background and hence 



there is no term corresponding to the first term in Eq. (2.24). In the appendix of Ref. [25], it has been shown that 
the expansion should be 

P..(?,n('^) - -^..(n('0^/BB (|.) + V(n«)P,^/BB (^) , (2.25) 

which is just a consequence of the fact that there is no background polarisation. We will check that Y"^^ vanishes at 
first order as it is not generated by collisions at this order. Similarly to the expansion of the intensity part, we want 
to expand the distribution function around a Planck spectrum at the background temperature T rather than the local 



number density temperature T. Thus we expand Eq. (2.25) in <d up to first order to get 



P^. = -(1 + 3e)7'^.^/BB (^) + (V + eV^.)VpBB (^) . (2.26) 

Again here, in order to simplify the notation, it is implied that Q, 7^^^ and Y"^^ depend on x^ and n^^'. 

III. BOLTZMANN EQUATION 

Now that we have all the tools at hand, we are ready to formulate the Boltzmann equation for polarized radiation 
in the cosmological context, and extract its spectral components. This section is entirely dedicated to this task. Given 
that the complete and detailed derivation can be rather lengthy, all details which are not necessary in a first reading 
are gathered in Appendix [Cj We first present the general expression of the Boltzmann equation for a tensor- valued 
distribution function. Since the Boltzmann equation is the description of how this distribution function evolves along 
a photon geodesic, it is necessary to perturb the geodesic equation up to second order. We then show how the 
Boltzmann equation can be split into its main spectral components, that is into a temperature and a distortion. 
Finally we write the explicit forms of the free-streaming part and the collision part of the Boltzmann equation. 

A. Boltzmann equation 

The evolution of the tensor-valued distribution function is dictated by the Boltzmann equation [28j 



where V/VX is the covariant derivative along a photon trajectory x''(A) and the momentum and C^j^ is the associated 
colhsion term. The exphcit form oi'Df^^/'DX is 



P/p. _ „ . dxP 5i^dg« 
PA - "•'^'^ dA ^ 9q(^) dA ' 



(3.2) 



where V^ indicates a covariant derivative associated with g^^. Using spherical coordinates (q, n*^*)) instead of g*^'-* for 
the momentum space, the Liouville operator, that is the l.h.s of Eq. (3.1 1, reads 



C PC rr 



VX 



^n i>v ^ Tlpa 



daf^ 9/^,. ding 
dA d In q dA 
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dn(*) 



dA 



(3.3) 



where Drn is a covariant derivative for momentum. The detail of the construction of such derivative is discussed in 



'[i] 



Appendix C 2 



equation can 



Thanks to the operation of the projection S^J'Sy" onto Vfp^/'DX, the left hand side of the Boltzmann 
30 decomposed into the /, V and P^^ parts similarly to Eq. (2.181; 



S.fS,/^^ = \mS^,„ + L[V]p, + ;^L[y]ep^,,e(o)''7i° 



'' " VX 

where the corresponding Liouville operators are defined as 

VI 



m^ 
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L[V]^ 



vv 



VX ' ' '^-^ '^ *^ VX ' ' ' VX 

Here, the operator V/VX on a scalar distribution function /, takes the simpler form 



VI 
VX 



(x^g,nW) 



da;^ df ding dnW 



dA 9 In g dA 



dA 



In a similar manner, one can also decompose the collision term, that is the r.h.s of Eq. (3.1) as 



(^iiiy = xL- S^iy + C„j^ + -6 ep^tiyo-ero) n 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



so that after extracting the trace, symmetric traceless and antisymmetric parts we get obviously L[I] — C^, L[P]fj_^ ~ 
C^^ and L\V] = C^ . Beware that this does not mean that the intensity, linear polarization and circular polarization 
evolve independently, since for instance C^ is the " intensity part" of the collision term but it may involve in general 
all components /, P^j^ and V of the tensor- valued distribution function. As a matter of fact, Compton collision does 
indeed intermix intensity and linear polarization, whereas circular polarization evolves independently. 



B. Geodesic equation and momentum evolution 



From Eq (2. 12a I and the definition of momentum, dx^ /AX = p^ we obtain 



dA a 
dx' 
dA 



= 12(1^")^ 



4 (n(*) -P^- h^,n^^^) , 



(3.8) 
(3.9) 



where only the first order terms are kept. In terms of (g, «.*■*'), the geodesic equation leads to the evolution equation 
for the conformal energy 



ding 
dA 



)^,^n(*^ + /3,,jn(*)7i(^) - K/n^'^n'-^^ + afa,,^^*) - l3,^jn'^'^n^^^ + h./n^'^n'^^A 



As for the direction evolution, up to first order in perturbations, we obtain 



U) 



dA 



^(«)0-) L^^. - (/?fc^^. - /i,fc')n(^-) + ih,i,k - hkiM'^'n^'^ 



(3.10) 



(3.11) 



It is obvious that n(^i-)dn^''^ /dX = 1 as it ought to be since n^"^^ is a unit vector. We need these expressions only at 
first order, except for the evolution of q because the background distribution function is constant in space-time (see 
below). 

Before closing this subsection, we mention that wc are free to choose another afhnc parameter than A, to label a 
point on a geodesic. A convenient choice is to take the conformal time 77 at each point of space-time crossed by the 
geodesic. The advantage of such choice, is that for photons having the same direction (and which thus follow the 
same path), but not the same energy, the same conformal time 77 would correspond to the same point of the geodesic. 
We can trade A for r] using Eq (3.8 1, that is with 



dr/ 
dX 
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a^ arj q 

The evolution of position, conformal energy, and direction, take then the form 
dx^ dx^ dA 



di] dA di] 

d In g d In q dA 

drj dA di] 

dn(*) dA 



dn(*) 
drj 



dA dry 






(3.12) 

(3.13a) 

(3.13b) 
(3.13c) 



C. Spectral decomposition of the Boltzmann equation 

We are now in position of writing down explicitly the Boltzmann equation, expanding the orders of perturbations, 
and separating the spectral components. Let us first look at the formal structure of the Boltzmann equation, especially 
focusing on the spectral decomposition. At the background level, the Boltzmann equation yields 






dln(aT) d/BB(a;) 



dA 



dlni 



0. 



(3.14) 



x=q/(aT) 



This implies that I has no time dependence and T scales as 1/a. One can conclude that the Planck distribution does 
not change in time if the initial distribution is given by the Planck one. This does not mean that the radiation is 
not losing energy as the universe expands. Indeed, since the physical energy of a photon is not the conformal energy 
q but q/a, then p ex J I{q/a)^dq/a ex a~* as expected. This background result for the scaling of T is useful as it 
implies that only the partial derivative with respect to q on I-bb[q / {o,T)] are relevant, and this motivates our use of 
the conformal energy. 



Now the action of the Liouvillc operator on the intensity, Eq (2.24), is expanded up to the second order as 
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(3.15) 



where we used dlBB/drj = 0. 



When we compare this formulation of the Liouville operator with the spectral decomposition (2.241, we are tempted 



to say that the expression inside the first square brackets contributes to the evolution of the temperature due to free- 
streaming, and that the expression inside the second square brackets is very closely related to the evolution of the 
distortion. It order to give a clear meaning to this assertion wc decompose Eq. (3.15) according to 



L[I]^ 



q 



- c 



2ec'^] 



dh 



BB 



/ 9 In g 



/:^ + 9/:®)2?2/, 



BB 



(3.16) 



This spectral decomposition is motivated by the fact that i) in the case where the conformal energy is not affected by 
free-streaming, dlnq/dA = 0, then {q/a^)C^ simply reduces to L[9]; and ii) the prefactor q/a^ is introduced because 
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the Liouville term is expected to be proportional to g/a^, as it can be inferred from the expUcit form (3.81 of drj/dX. 
From a comparison of Eq (3.151 with this decomposition, we have 



'Lc^ = m-ii + e)^^ 



a 



dA ' 



m 



(3.17) 
(3.18) 



We must bear in mind that in these expressions, even though the operator L[.] has been defined in Eq. (3.51 for 
functions of {ri,x'^,q,n^'^'), it is applied on the spectral components Q and y which do not depend on q. 

Since the Liouville operator is equated to the collision term in the Boltzmann equation, it is convenient to decompose 
the collision term in the same manner as the Liouville term. That is, it is decomposed as 



C 



dh 



(c« + 2ec«) ^ + {c^ + ec'-'ypBB 



(3.19) 



such that the spectral components of the Boltzmann equation can be formally very simple and are given by 



(3.20) 



This decomposition means that once the spectral decomposition of the collision term is known (C and C^), then we 
only need to obtain the spectral decomposition of the Liouville term from Eqs (3.171. 



We follow the same logic for polarization. First, the corresponding Liouville operator reads 

^■2 1 



L[P]^, = -L 



(1 + 36)7'^, 



9Ibb_j, dlnqcK/BB 
d\nq ^^ dA dlnq"^ 



Y., 



(l + 3e)i[7']^, + 3 L[e]- 



dA '' ^^ 



OIbb 



-DglBB 



dlnq 



L[Y] 



fiiy 



im -^w, 



^^ ,. ^.^ eL[PU 



T^qlBB ■ 

(3.21) 



For the same reasons as in the case of intensity, it appears natural to decompose this Liouville operator into spectral 
components according to 



L[^U 



dh 



(1 + m)cl, + ic'^v,,] -^ + [cl, + c^v,. + ec^^)vpBB 



which implies that the spectral components are given by 

The collision term must then follow the same type of decomposition, that is 



2 ^fj-u — ^['^ \^ll' ■ 



(3.22) 



(3.23) 



CL = 



(l + 3e)C + 3C®7'^, 



f|| + (C + c^X. + eC)p,^/BB 



(3.24) 



so that, again, the spectral components of the polarized part of the Boltzmann equation take the formally simple form 









(3.25) 



Again, this decomposition means that once the spectral decomposition of the collision term is known (C^^^ and Cj^), 
then we only need to obtain the spectral decomposition of the Liouville term from Eqs (3.23) bearing in mind that 
the Liouville operator L[.] applies on functions which do not depend on q, but only on {rj,x^ ,rv'^>). 



D. Temperature and spectral distortion of Liouville operators 



Now that the spectral separation of the Boltzmann equation is performed, it is time to expand the equations 
obtained in orders of perturbations. In the next two sections, we present such expansion for the temperature and 
spectral distortion parts of the Boltzmann equation. The case of polarization is reported in Appendix |C 4[ 
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At first order in perturbation, with Eqs. (3.8) and (3.10), the Boltzmann equation leads to 



(3.26) 



Note that there is absolutely no q-dependence, nor scale factor a in this expression, meaning that our spectral 
decomposition performed in Eq. (3.16) is adequate. Concerning the spectr al dist ortion part, the Liouville part at 
first order is C^ = y' + y,in^^\ but since C^ = at first order (see section HIE I, one can conclude that only the 
temperature part evolves at first order and no spectral distortion is induced. 

Up to the second order, the Boltzmann equation for the temperature is given by 



£« = 



L[e]-(i + e) 



ding 
dA 



e' + e,in(') + a,,n(*) - (3,.n^''U^^^ + hu'n^^Ji'^^^ 



a 

q 



drj 
dA 



(1) 



e' 



dx' 
dA 



(1) 



e. 



dn(») 



dA 



(1) 



D,e - 



ding 



dA 



(l)x(l) 



-e 



ding 



dA 



(!)■ 



= c' 



As for the spectral distortion, the Boltzmann equation is given, even at second order by 



(3.27) 



(3.28) 



This means simply that gravitational effects do not induce spectral distortions, and this result holds actually non- 
perturbatively. 

Before ending this subsection, for the sake of completeness, we shall write down the most general form of the 
second order Boltzmann equation for the intensity. Using Eqs. (3.8), (3.10), (3.9) and (3.11), the detailed form of the 
evolution equation for temperature obtained in Eq. (3.27) is 



- aQ' - {(3' + /i',n(^))e,, - [a,, - {13^., ~ h,/)n^^'^ + (h^j - h.^M'^^^''^ 

- (aas + a^jhii'^n'^'^ - a(-/3,j + hi/) + (3''h,j,k + (/3fe,i - ft,fc + 2h^k)h^ 



D*e 



rSi)nU) 



(3.29) 



E. Temperature and spectral distortion of collision terms 

The expression of the collision term has been derived by taking into account only the intensity in [12l [29] and then 
it was extended to include the effect of polarisation in [T31 [HI [TT] . Here we summarize the result obtained by Benekc 
et al. jl7] applying the decomposition of the distribution function into intensity and linear polarisation. 

The complete expression of the collision term for intensity is given by 



C® ^an.ar ( -6 + (6) - -5«(^) [(0"i(»)a)) - 2(7'(,)(,)) 
C^ ^an^ar (-y + (y) - \s^^^^^^ [{ym(^)U)) - 2(%(.)) 



+ wWnai + 



W-^5^ + 5«WQ^)(,.)+5.C^ 



-5- + 5«0-)2^^^^.^ 



where 5"^, 5^, Qj^-.,,-. and QX-ir,) s-i'G quadratic contributions defined by 



(3.30) 
(3.31) 



S' 



e^ + (e^) - 26(6) - 6i;«n(,) + 2(6)i;«n(,) - 2(6n(,))«« - -i;«z;(,) + (i;«n(,))2 , 



5^ = ^©' + \(^') - ©(Q) - Q^^'^^W + (0)"^*^"(^) - (en(,o)^'(*^ + ^«(^'«(.) + liv^'^n^:)^ , 



(3.32) 
(3.33) 
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^(i)U) 



= --(e2m(,)(j)) 



e(em(,)(j)) + - (e7i(j))f(i) + {en^i))v^j) 



,ik) 



%-)(Q^ 



Ht)(k)) 



'\^) 



(Or 



\j)(k)i 



- ^„(fe) 



2n(fe)(eTO(j)(j)) + (em(,)(j)n(fc)) 



- 5"(0«0-) 



+ o<Q^wa))-30(n.)a)) + o" 



(fc) 



,(fc) 



^(i){j) 



= -^(0^ 



(^(i)(fc)"-(j)) + {'P(3)(k)n(v, 

-e(eTO(,)(j)) - -t 

3 , 3 



3 



,(*) 



('^W0)'^(*=)) ^2n(fc)('P(i)(j)) 



i(0{i) 



n(fe)(em(j)Q-)) - (em(i)Q-)n(fe)) 



20 



V(€iV(j) 



+ ^(0^00-)) - i0(nOO-)) + i^'^'^ [^(fe)(n^)0)> - (n»)0)"(^)) 



and m 



{i){j) = n(i)n{j) 



(3.34) 



(3.35) 



(5(i)(j)/3. Note that we have introduced the notation [Q) — J^Q = Jd^n'^^^Q, which 
corresponds to a muhipole extraction that we do not perform exphcitly here. Note also that 6e = Sne/fie is the 
fractional perturbation of the baryons number density, and w'*' are the tetrad components of the baryons spatial 
velocity. Again, we also defer the expression of the collision term for polarisation to Appendix |C 4| 



IV. GAUGE DEPENDENCE OF THE DISTRIBUTION FUNCTION 



Now that we have established the Boltzmann equation, up to second order, with its spectral components separated, 
we investigate the gauge dependence of its constituents. Eventually the Boltzmann equation itself should be gauge- 
invariant, so if we are able to check explicitly that the Boltzmann equation is gauge-invariant, this means that it 
is very likely that i) the perturbative expansion of the equation is correct; and ii) the gauge transformation rules 
for all its constituents (metric and distribution function perturbations) are correctly understood. We thus consider 
this verification as a consistency test. This section is dedicated entirely to this task. We first review the gauge 
dependence for tensors, and deduce how it can be extended to a scalar distribution function. The case of a tensor- 
valued distribution function, even though it is the less trivial part, is treated in appendix |D 1[ We then infer what 
should be the transformation rule of the Liouville and collision operators, and in order to complete the consistency 
test, we check that the perturbed expressions of the Liouville and collision operators do indeed transform following 
these rules. 



A. Coordinates on the manifold 

We need to specify how the functional dependence of the quantities appearing in the Boltzmann equation (and 
in the Einstein equation) is obtained. If we consider a scalar function / : A^ i— >■ M on the space-time manifold A4, 
then the choice of a coordinates system ^ c : M"* i— )■ Al does not affect the geometrical meaning of this function, but 
it affects its functional form / o c : K'' i— >■ K, where o designates the composition rule, in the sense that for another 
coordinates system c : K^ i— J' Al , then / o c ^ foe. The confusion only arises from the fact that we often refer to / o c 
as / only. A distribution function / (that we take as a scalar-valued for simplicity here) is a function on the tangent 
bundle TA4 of the manifold and can be regarded as a function on manifold A^ describing the space-time and on the 
tangent space M"' (more precisely a restriction to the mass shell M.^) of each point. It is thus a function 



/ : (TM) ^ : 



(4.1) 



Again the particular choice of coordinates on the tangent space does not affect the geometrical meaning of the function 
but its functional form. Once a choice c of coordinates on the manifold A4 has been made, there is a natural basis, 
called canonical basis, that is made of the partial derivatives with respect to the coordinates. This leads to a natural 
coordinates system Tc for the tangent space at each point. Thus {c,Tc) : M'' x M'^ i— >• TA^, is a coordinates system 
for the tangent bundle. In order to simplify the notation we will note (c, Tc) as simply as c. 

Furthermore, in this paper we use coordinates in the tangent space described by the tetrad components. More 
specifically we use the components of the conformal momentum in spherical coordinates in this tetrad basis, {q,n^'^'). 



We assume for the simplicity of the argument that one system of coordinates is enough to cover the entire manifold. 
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There is a problem with such a choice since the tetrad basis is not unique. However, once a choice c of coordinates on 
the manifold A4 has been made, the tetrad basis might be completely fixed from the metric through a prescription 
described in § |II A| Once a coordinates system c has been chosen, and the tetrad is fixed thanks to this choice, we 
obtain the functional form of / in the form / o c : !&■* x M^ H' M. 

B. Geometrical interpretation of the gauge 

When performing perturbations around a background FLRW space-time, we need to have a one-to-one correspon- 
dence between the background space-time A4 and the physical (and perturbed) space-time A4. This can be completely 
defined geometrically [30 but we take a shorter approach. If we have two sets of coordinates ^ c : M'' i— )• 7W and 
c : M^ i~> A^ on the background and the perturbed space-time, then we identify points with the same coordinates, 
that is, we identify points with c o c~^ : M. i— >■ M.. Since we could have chosen different sets of coordinates, there is 
some freedom in this choice, which is known as the gauge freedom. 

On the background, the symmetries can justify that we can find a preferred choice of coordinates. For instance 
for a flat FLRW space-time within a given background cosmology, it is enough to choose that the time coordinate 
is the proper time of observers with 4-velocity orthogonal to the homogeneous surfaces, that is the proper time of 
comoving observers. On a given homogeneous surface, there are preferred choices of Cartesian coordinates, since it is 
conformally related to E^, and all these Cartesian systems on the spatial homogeneous surfaces are related by global 
translation and rotation in M? which are irrelevant given the homogeneity. And once a coordinate system has been 
chosen on a homogeneous surface it can be Lie dragged by the comoving observers to any homogeneous surface. So 
essentially there is a unique mapping c from M^ to the background manifold M. This point is illustrated in the left 
part of Fig. [T] 

However on the physical space-time we could consider another coordinates system c : M"* n- A^, and this leads to 
a different identification through co c~^. The fact that we fix the system of coordinates on the background space- 
time but there is still some freedom on the physical space-time for the choice of coordinates leads to a freedom in 
the identification between points of these two space-times. With c, a given point P G M would be labelled by the 
coordinates x^, that is c{x) = P, and with c it would be labeled by the coordinates x'^, that is c{x) — P (see the 
upper part of Fig. nj). For every point, there exist four numbers ^'^ = (T, L') such that 

S:^'{x) = x^' +C{x). (4.2) 

We should note that, although there is an index in the notation, <^'^ is not a vector field on Ai nor on M, but it can be 
seen as a vector field on W^. In the literature, a vector field (" is often used to generate a coordinate transformation 



?;m 



1 



X 



exp (£c) x" = a:'' + C" + t^C" A" + ■■■ ■ (4.3) 



2 



In this paper we adopt the former definition Eq. (4.2 1 and care must be taken when comparing our transformation 



rules with those in the literature. Note that in the rest of this section, we will use extensively the notation 

x = x^, X = x^' , (4.4) 

even though a;'' and x'^ are not vectors but just coordinates. 

C. Metric transformation 

In general the coordinate transformation rule of any tensorial quantity is given by 

Tod{x)=Toc{x). (4.5) 

This means that it is invariant under a coordinates transformation since it is geometrically defined and independent 
of the coordinates used to parametrise the manifold. From this we can deduce the gauge transformation of its 



Again here for simplicity, we assume that such coordinates system covers the whole manifold. 
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X 



FIG. 1: In this figure, we represent all four dimensional spacetimes with only two dimensions. First, we have noted that there 
is essentially a unique way to map R'* to the background manifold, that is to relate the top-left to the bottom left. However, 
there are several ways to relate R'* to the perturbed manifold, and consequently to relate the background manifold to the 
perturbed manifold. We represented a coordinates system c and another coordinates system c which relate R** (top-left) to the 
physical manifold (top-right). For each point P of the physical manifold, there is a set of four numbers x'^ and another set of 
four numbers x'^ such that c{x'^) — c{x^). Furthermore, each coordinates system has different surfaces of constant time, and 
thus different set of tetrads, given that the null tetrad is always chosen to be orthogonal to the constant-time surfaces. When 
using the tetrad field to extract the components of a given momentum (that is a point in the tangent bundle), this will lead 
to different components, depending on the coordinates system chosen, but since this is the same momentum at the same point 
and the tetrad are normalized, we can relate the components by a Lorentz transformation. For a given point of the tangent 
bundle, that is, for a given point of the manifold and a given momentum, we always have c{x^ ,p^^') = c(i^,p'"'). 



components, which is defined as the transformation when the tensors are compared at the same coordinate in two 
different coordinates systems. Let us consider in particular the metric. Since g^^ = g{d /dx^ ,d/dx'^) and gfii, = 
g{d/dx'^, d/dx'^), we then deduce the usual coordinate transformation rule of a 2-form 



.g/ii> o £{x) 



dx" dx^ 
dxt^ dx'^ 



gap o c{x) . 



(4.6) 



Expanding the left hand side which is evaluated at x around x leads the gauge transformation rule up to the second 
order 



1^2, ^ ^ 1 



9&P ° c(a;) = gap o c{x) - C^g^p o c{x) + -C^g^p o c{x) + -C^^g^p o c{x) 



(4.7) 
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where ^^ designates the quantity £,'^,v£,'^ ■ We emphasize that the coordinate transformation (4.6 1 is the transformation 



of components at the same point of space-time, whereas the gauge transformation (4.7) is the transformation of the 



components at different points which have the same coordinates in two different coordinates systems. Since the 
notation can become rather cumbersome if we specify which coordinates system c or c is to be used, we shah indicate 
it only when it is the new coordinates system c. 

Throughout this paper, we use the symmetrization and anti-symmetrization definitions 



^{^.f) = 2 ^^^^'' ^ ^''t^^ 



X 



[pH ~ 9 ^ f^'^ ^^m) 



(4.8) 



Up to second order, from the transformation (4.7) and the decomposition (2.1), we obtain the gauge transformation 
for the ADM variables as 

a^a-HT-T' + -{H^ +n')T^ +n{2T'T + T,L') - HaT 



aT' - a'T - a^,V + /S'T, + T"T + T' 



rrv 



-T-'T, 



pt ^ P^ j^ rp,i _ j^i' ^ 2aT''' 
- (2T'T^* + TT'^') ^ 



2hij = 2h,j - 2'HT5,j - 2L(,j) - iTiTh^j 






TV" + T'V + TjL'-J 

[2%^ + 'H')T'^5, 



(i\jV - 2h'^Tj 



T'\jV 



I/yV 



2H(Tr + TkL'')S,_ 



-^/^(i^j) ~ 2/iij T — 2hij^kL — Ah(i\kL _|j) 



TaT^j + 2TL(^ij-)' + 2T(^,jij)' + 2L(.jj)feL ' + 2L(j|^fcL 'jj) + LkAL j 



(4.9a) 
(4.9b) 

(4.9c) 



These relations must be understood as follows, a^a at first order is the 00 component of the first order metric in 
the c coordinates system, and taken at the point of coordinates x, and is thus equal to — Iffop ° c{x). Instead, a^a 
which means a'^d{x), when considered at first order is the 00 component of the first order metric in the c coordinates 
system, but also taken at the point of coordinates x, that is n'-^nn 



Ifffifi^ ° c(a;). 



D. Tangent space basis and tetrads 



As discussed in Section II C the transformation of the basis on the tangent space is entirely linked to the coordinates 
change in the base manifold. Usually, the canonical basis d/dx^ and the corresponding forms da;^ are used as a basis 
of the tangent space and they transform according to 



dx^ d 



dx'^ 
dx'^ o c(x) = - — dx^cix) . 



(4.10) 



When we consider the components of the metric, these components refer to this canonical basis. However we will 
use the tetrad field as a basis for the tangent space. We thus need to relate ej^) o c{x) with e(a) ° c{x) in a similar 
fashion. Since this is a relation between two orthonormal basis at the same point of space-time, there exists a Lorentz 
transformation ^ L\ such that we can relate the two tetrad fields associated with the coordinates systems c and c 
(see the bottom right part of Fig. IT] for an illustration of this) as 



;;^^^ - A i'')/ 



£(") nn(^\ - A^'') 



,A(^) A^'') -• 



e(a) o c{x) = A(^^) >{x)e(^h)ix) , e^"' o d{x) = A^ lj,)ix)e^ '(x) , ?7(c)(d)A^ '(^^A^ ^^^ = 7]^a){b} 

The components of this Lorentz transformation are given, up to first order, by 

(0) ' (i) (0) * ' 0) 0) .j] 



(4.11) 



(4.12) 



At second order, it proves more useful to relate the components of these quantities, that is, to relate e(a)^ o c{x) = 
dx'^ o c{x)[e(^a)] to e(Q)'^(a;) = dx'^{x)[e(^a)]- For the former components, we must use (2.6) with a(^), a{x), $i[x) and 



hij{x) [that isgnp o c(£)] which can be deduced from the rule (4.9) just by shifting the argument of the left hand side 



of the rules (4.9) from x to x. For the latter we must use (2.6| with a{rf)^ ol{x), /3i{x) and hij{x) [that is 5^j,(a;)] 
We do not report the corresponding expression since we will work instead directly with the perturbation components 
of the metric in the next section. 
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E. Momentum, energy q and direction n 



(i) 



The components of the momentum of a particle in the canonical basis transform as 



p^ o c{x) 



dx" 



p^{x). 



with p''EEdi^(p), p^ = dx^(p). 



(4.13) 



However, we are going to use the tetrad basis in the tangent space rather than the canonical basis. We thus want to 
express p^°'^ o c{x) = e^"' o c(ai)[p], [or qo c{x) ^ and n^*-* o c{x)] as a function of ^("'(0;) = e'--°-Mx)\ n\ \ot q(x) and 
n^*^(a;)]. Eventually, we will prefer to use q and n^*^ rather than p^'^K From the definition of q, Eq (2.13a) we then 
obtain the desired transformation relation 



q o c{x) = q{x) + Sq{x) = q{x)[l + S In q{x)] , 



(4.14) 



q o c{x) = a {f))Np o c{x) 



= q{x) 



l+'HT + T,n^'^ + -{'H' + 'H^)T^ + HTT^n^''^ ~ T'T,n^''> + -T^'T^ + ^,» f"^^*^ - h'^n'^^A 



As for n^*' from Eq (2.13bl, its transformation 

nF> o dix) = n^''> (x) + Sn^''> (x) 
is given by 



.W 



c{x) 



P' 



1 

N 



W 



PL 

" pO 



,w 



(x) + S'^'^^^^Tj + L^''^^1 



c{x) 



ij) ■ 



(4.15) 



(4.16) 



(4.17) 



Here when the argument is not specified, it is (a;). Note that for future use, we have defined in these expressions the 
differences 6q, ^Ing, and Sn^"^' . The first order and second order perturbation of these can be read directly from the 

expressions above. We also define Sq^"^^ = q {S\nq)n^'^^ + Sn''^^ . 

It is worth stressing that these differences are measured at the same point. For instance (5n*^*^(a;) — n^*) o c{x) — 
n'^^^x), and they do not vanish because in one case we use the tetrads e^*) associated with the coordinate system c to 
obtain the components, and in another case we use the tetrads e^*' associated with the coordinate system c. It is the 
basis at a given point of space-time that changes when we change the coordinate system, not the momentum itself. 
This point is illustrated in the bottom right part of Fig. (II]) 



Scalar distribution function 



If we were using the natural basis associated with a coordinate system (the canonical basis) for the tangent space, 
then any scalar function on the tangent bundle TM , that is a function of the space-time position and of the tangent 
space at each point, would transform as 



Io£{x,p^') =I{x,p^'), 



(4.18) 



where x and x are related by (4.2) and p^ and p'^ are related by (4.13). Again this rule is a statement that the 



function is invariant under a change of coordinates because it is defined purely geometrically. 

However, as mentioned earlier, we use the basis of the tetrad field, e(j) and e^*^ to obtain the components of 
momentum not the canonical basis. The tetrads are also completely determined by the choice of coordinates due to 



Under our conventions, q = ap^-^' . 
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our prescription (|2.6|. The tetrad field, though being of tensorial nature, is not invariant as in Eq. (4.5 1. Furthermore, 



as mentioned earher, we also work with the conformal momentum q rather than the momentum itself p. Given this 
choice for the basis of the tangent space, the scalar function transforms as 



/oS(£,gW)=/(a;,gW)^ 



Wl 



(4.19) 



that is, it is unchanged when it is evaluated at the same point of the tangent bundle. On the other hand, the gauge 
transformation is a transformation rule at the same coordinate point and it is the relation between / o c(a;, g*-*') and 

I{x. g^'-*). We then need the expressions of g*^*' o c(x) in terms of q'^^\x) in spherical coordinates, which are derived in 
the previous section in (4.15). At first order, using the fact that the background distribution function cannot depend 
on the direction n*^*^ , we obtain 



d 



I o c{x, 9, nW) + S^t'd^ + Sq— I o c{x, q, n^) = I{x, q, n'-'>) 



,W^ 



,W^ 



(4.20) 



Given that the background distribution function also depends neither on time nor on space but only on q, we obtain 

dliq) 



/(i)o2(a?,q,nW) = /(i)(a;,g,n('))-(51n( 



91ng 



At second order, we obtain 

/ o £{x, q, n«) + (ff^d^ + Sq-^ + 5n^^D(,) + \h^ ^\ I o ~c{x, q, n^) - I{x, q, n^) . 
Using the first order expressions, the gauge transformation rule reads 



(4.21) 



(4.22) 



/(2)o5(a;,q,n(*)) 



ll''\^, 



q,n 



W^ 



^^'d.+Slnq 



d\nq 



+ (5n(*)D(,))/(i)(a;,g,n(^)) 



{^'9^ + ^-'''9^)(^'-'i^d^,^('i^ + l(^'-'i^'(^l-^ 



dlnq 



I{q) , (4.23) 



where we used dSlnq/d\nq = 0. 

This transformation rule can be applied to / or V^ since these are scalar valued distribution functions. However we 
are interested in the transformation rule for the spectral components of /. Using the decomposition Eq (2.24) for / 



we obtain that the temperature is transforming under a gauge transformation as (noting for simplicity Q = Q o c) 

e(i)^e(i) + (51nq)(i), (4.24) 

le(2) = le(2) + i(51n<7)(2) + e(i)(51ng)(i) - f ^f^ ^ + Sn^''> D ^,^ \e^^'> + {Slnq^^A , (4.25) 

where it is implied that all quantities are evaluated either at x or at (a;, g, rv '>\ The detailed form of the transformation 
rule can then be obtained just by considering the perturbations of q and n*, (5 In g and 8'rv-'^' , which have been obtained 



in Eqs. (4.15). For completeness we report it here 



e = e + -HT + T,n^'^ + -{n^ - 'H')T^ + nT{T,n^''> - T') - {T'T, + TT'^)n^'^ 



+ {HT + T^,n'^'^)Q - e^e,^ - (s'(')(^)Tj + L^\J^^n'^^^\ D (^^ . (4.26) 

Finally, the gauge transformation of y is trivial. Since y vanishes at first order, y is gauge invariant at second order, 
and it can also be checked directly by extracting y out of the transformation rule of / at second order (4.19). 



G. Baryons fluid description 

In order to obtain the complete gauge transformation of the collision term, we need the gauge transformation rule 
of the baryons fluid velocity in the tetrad frame up to second order, since it appears in the collision term. We obtain 



z;^'^ o c 



(a;) = v^'\x) + T'' + aT' - h'jT-^ - T'T'' + L^\j]iv^^'^ + T') - T{v^'^' + T^'') - V (y^'^ + T'^,, 



(4.27) 
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We also need the gauge transformation rule up to first order of the electrons density and it is easily obtained to be 

Seo£{x)^S^{x)+ iHTn^ . (4.28) 

H. Gauge dependence of the Boltzmann equation and Gauge invariant form 

Having derived all the necessary gauge transformation rules, it is now possible to check the gauge dependence of 
the derived second order Boltzmann equation and collision term explicitly. More precisely we shall check that the 
Liouville and the collision terms of the Boltzmann equation transform as they should do. The Liouville term and 
the Collision term are distribution functions (scalar or tensor valued depending whether or not we are considering 



the intensity or pola risation). From the transformation rule of a scalar distribution function (4.23) and the spectral 



decomposition (3.16) we deduce that C^ and C^ are gauge invariant and £® and C® shoulc 
second order as (noting £® = £® o c and C^ = C^ o c 



transforms up to the 



£6 ^ £0 _ /^M^^ ^ 5n(»)D(,,) V® (1) _^ 2HTC^^ (^^ , (4.29) 

Ce = C® - (ed^ + <5n(')i:'(,))c®(i' + 2HTC'^^^^ . (4.30) 

After very long and tedious but straightforward calculations using all the transformation rules derived so far for 
the distribution function, the metric components and the baryons velocity and energy density, we have checked that 
the Liouville operator C and the collision term C® actually transform as in the above equations. This completes 
the consistency test of the Boltzmann equation that we have derived as well as all the gauge transformation rules 
obtained for its constituents. 

The same property is found of course for the circular polarisation since in that case the collision term vanishes. As 
for polarisation, we have also checked that the Liouville operator 'C^w ■■, and the collision term C^^w ■, transform like 
tensor valued quantities (see the details in Appendix |D2[) . More importantly, we have checked that the Liouville and 
collision terms for the spectral distortions (£/^w ^ andCX^w J are gauge invariant as it should be since they vanish on 
the background and first-order spacetimes. 

The gauge invariance of the Boltzmann equation, as in the case of Einstein equation and in general for covariant 
equations, enables us to write it down in terms of gauge invariant variables. In practice, it is equivalent to completely 
fix the gauge and write down the equations in term of the perturbation in this gauge. 

V. SUMMARY AND DISCUSSION 

In this paper we derived the second order Boltzmann equation in the most general manner incorporating polarisation 
and without fixing a gauge. In order to describe the polarisation of photon, we used a formalism based on a tensor- 
valued distribution function. We performed the separation between temperature and spectral distortion for the 
intensity and we also extended this separation to polarisation. 

We then derived the gauge transformation rules for the metric, the momentum and the distribution function to see 
how those quantities are mixed under the gauge transformation. As an application, we checked the gauge dependence 
of the derived Boltzmann equation under a gauge transformation and obtained consistent transformation rules. This 
is a non-trivial check of the correctness of the derived equations as well as the gauge transformation rules. 

We now discuss two issues related to the gauge dependence in the Boltzmann equation. 

A. Gauge dependence of lensing term 

It is well known that the lensing term 

dn^') 
£ D DiQ : (lensing term) , (5.1) 

dA 

which is written in terms of the conventional lensing potential in the Newtonian gauge significantly affects the bis- 
pectrum of CMB [2l[3Tl[32]- Indeed, the correlation between the lensing and ISW effect is the dominant contribution 
to the bias for the local-type non-Gaussianity in Planck. However it is very hard to include this contribution in the 
line of sight integration and evaluate it until today. Usually, the lensing effect is added separately to the final result 
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obtained in the Poisson gauge. However the effect from this lensing term depends on the gauge choice. Actually, as 
we have seen above, the lensing term is mixed with other terms under the gauge transformation. This means that 
some of lensing effects in a specific gauge are absorbed into other effects in another gauge. In principle, there exists 
a gauge where we can avoid the difficult computation of this lensing term to some extent by evaluating other more 
tractable terms. Since we have derived the Boltzmann equation without choosing any specific gauge, it should be 
possible to investigate this possibility further. 

B. Observed temperature anisotropies 

Here we make a comment on the observed temperature anisotropies. In the main part of this paper, we have shown 
that the second order Boltzmann equation is gauge invariant and thus it can be written in terms of gauge invariant 
quantities. However, there is a subtlety in the meaning of "gauge invariance". This originates from our choice of the 
local inertial frame. 



As is clear from Eq (2.6), we always choose the local inertial frame so that the three-velocity vanishes, -0* = and 
there is no rotation of the spatial axis relative to the background spatial coordinate axis, let us call Oi = 0. In order 
to achieve this, the local inertial frame has to be changed when we perform a gauge transformation. If we were to 
identify this local inertial frame as the one of an observer, we would be lead to consider different observers in different 
gauges. This is clear from the gauge transformation at the first order: 

e ^ e + -HT + T,ri\ (5.2) 

The last term comes from a change of local inertial frame. By fixing the gauge we can promote © to the gauge 
invariant temperature fluctuations but these are temperature fluctuations observed by an observer with w* = and 
9i = in this gauge. In order to evaluate temperature fluctuations observed by a different observer, we need to change 
the local inertial frame. Alternatively, we can perform a gauge transformation keeping the conditions w' = and 
9i = for the local inertial frame so that this frame coincides with the one of the observer. 

In all the literature, the second order temperature anisotropies are calculated in the Poisson gauge so far, with a 
specific choice of the local inertial frame. Strictly speaking this is not the temperature anisotropies that we observe 
as there is no reason for us to be comoving with the local inertial frame associated with such a gauge. One thus needs 



to change the local inertial frame or change a gauge. At first order, this was not an issue. As is clear from (5.2 1, 
the change of gauge and the local inertial frame only affect the monopole i = and dipole ^ = 1 if we expand 
the temperature anisotropies into multipole components. Thus, the £ > 2 modes are not affected by the change of 
observers. However this is no longer the case at the second order. In the second order gauge transformation, there 
are terms that are convolutions of the first order temperature anisotropies and the gauge transformation; 

e^&-Ce.^,-Sn(''>D(^,-,e + --- . (5.3) 

These terms affect the observed temperatures even for the £ > 2 modes. 

In order to define the " observed temperature anisotropies" , we should keep the conditions v^ — and 9i = for the 
local inertial frame and specify the gauge so that this local inertial frame coincides with our local inertial frame where 
we perform experiments. Thus special care must be taken when we compare theoretical predictions to observations. 
Our formula for the gauge transformation will be useful to investigate this issue further. 
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Appendix A: ADM variables and usual perturbation variables 



From the form of the metric in the ADM parametrisation (2.1 ), and the pertm-bation of the lapse function and the 
spatial metric given in the equations (2.2), the perturbed metric is expressed up to second order as 

ds^ = a^{T]) [-(1 + 2a + a^ - A/3')d?7^ + 2(% + 2hij)P^dx'dr] + (J^ + 2h,j)dx'dx^] . (Al) 

This has to be compared with the usual parametrisation of the perturbations of the metric which is in the form 

ds^ = a^ir^) [-(1 + 2A)djf + 2Bidx'dTj + (% + 2C^j)dx'dx^] , (A2) 

where C'ij can be further split into 2 scalar, 2 vector and 2 tensor degrees of freedom. By a direct comparison of these 
two parametrisations, the relation between the two metric parametrisations is 



500 : 1 + 2a + a' - ft/3* = 1 + 2A, 
9ij ■ "ij ")" 2hij = Oij + 2Gij . 



At first order we obtain 

A(i) = ad) , 
and the two parametrisations are the same. However, at second order we get the relations 



B^^-P^, 



rW _ 1,(1) 



A(2) = a(2) + a(i)2 - /3f )/3(i)* , Sf ) = /3f ) + 4/.(J)/3(i)^- , c'lf = h^ 



(A3a) 
(A3b) 
(A3c) 



(A4) 



(A5) 



Appendix B: Construction of the distribution function for polarised light 

We consider a two-dimensional polarisation plane defined by two unit complex vectors im and e(ii), which are 
mutually orthogonal, e(A)*^e(B)'^5A"' ~ '^(A)(B)- Any polarisation e can be represented by a superposition of e(^) in 
the form 



E 

A=I,II 



e^e(A) 



(Bl) 



The orthogonal vectors e(^) define a polarisation plane and we choose them to be orthogonal to the direction of the 
photon n^') and to the observer velocity ejo), 



HAfn g^^ = e(A)'^e(o) g^,^ 



(B2) 



We can also associate canonically polarisation forms through e' '^ = g^i/e-iA)'^i E^nd they will be also complex unit 
forms and mutually orthogonal. The polarisation density matrix /ab is defined so that the expected number of photon 
in a phase-space element with a polarisation state e is given by 



*A,B 



/(x,p,e) = /As(x,p)e* e 



(B3) 



With such a parametrisation, all the electromagnetic gauge degrees of freedom have been fixed and we parametrise 
the physical degrees of freedom of this density matrix by the usual Stokes parameters as 



/ab 



\(I + Q U-iV 
2[u + iV I-Q 



(B4) 



where it is implied that /ab and the Stokes parameters depend on the position x^ and on the momentum p^ [or 
(g, n'*-*) in spherical coordinates]. Jab is a Hermitian matrix since Jab = Iba- 

From the four-dimensional point of view, the polarisation density matrix is a tensor-valued distribution function. 
It is a 2-form defined by 



ftiK = I ABC* 



{A) (B) 



(B5) 
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and the expected number of photon in a phase-space element for a polarisation state e is given by 

/(x,p,e) = /^,(x,p)6*'^6''. 



(B6) 



This can be viewed as a multipolar expansion in the polarisation state. From (B2|, the tensor- valued distribution 
function is a projected quantity such that 



J^v "^/^ ^fi Ja(i 



(B7) 



It is then straightforward to realize that it can be decomposed according to (2.18) 



Appendix C: Boltzmann equation for the tensor- valued distribution function 

1. From a scalar valued to a tensor-valued distribution function 

In this section, we explain in detail how the Boltzmann equation for the tensor-valued distribution function can 
be obtained from the Boltzmann equation of a scalar distribution function. Since this scalar distribution function / 
depends on x^jp^ but also on e^, the action of the Liouville operator is given by 



V Ax" df dp" df de" df 

PA^^'''''^''''"^ = —^^ + — ^ + — ^ = C[f] . 

In the geometric optics approximation, p^ and e^ are parallel transported and we obtain 

X>p" Ap" 










VX AX 
^-vx- AX +h',^P ■ 



Using ( B6 ) , the term involving the evolution of polarisation is obtained as 
de" df 



d6" « Ae*f> 






dA 9e" " ^"'^ dA ^ ' ^"'^ dA ^ " .^5t^.,a/^c x^^c 
Combining this result with the space-time derivative term of the Liouville operator, we get 



P 



dx'-^ AXde'- "^ lV^/„^je 



and thus the Boltzmann equation (CI I can be rewritten as 

Vf 



- <:*' r." 



VX 



e" KV„ L, 



AX dp°' 



*•' = C[f] = e^C^.e*-- , 



(CI) 

(C2) 
(C3) 

(C4) 

(C5) 

(C6) 



The last equality defines the tensor-valued collision term. If we do fix the electromagnetic gauge condition for the 
collision term in the same manner as what we did for /^i^, that is, if C^p — S^/^Su Caji, then the Boltzmann equation 
for /^^ is given by 

Note that the use of the projectors is required because the components of the equation which are not in the polarisation 



plane are not fixed by (C6) 



From the canonical basis to the tetrad basis 



In the equation (C7), the Greek indices refer to a given coordinate system and its canonical basis for the tangent 



space, and the distribution function is a function of {x^,p'^). If we want to use instead an orthonormal basis for the 
tangent space, that is, to use the components p'*' = e^^'^p^ or the conformal momentum components q^^' = ap^^' in 
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the tetrad basis, then the Boltzmann equation can be modified accordingly. In order to do so, we need to be exphcit 
about the partial derivatives to emphasize which variables are to be kept constant when the partial derivatives are 
evaluated. The Boltzmann equation reads indeed 



Q an 13 T^faP _ c a q P I ^-y djafi 



vx 



dxi 



p''^Jsp-p''^pfo.5 + 



-yyS 



dX dp^ 



(C8) 



Using the properties 



dfa0 



dxi' 



dfap 



dxf^ 



we then deduce that 






dfal: 



dq(i) 



ae 



(i) 



dq 



(«) 



dA 



ae 



W "F „M ^("^^'V) 



dp^' 



dX 



dX 



(C9) 



q aq P 



VX 



S"S,^ p'< 



5/< 



aP 



dx'y 



-p-'T'^Js0-p''r'^pUs + 






dg(^) a/< 



aP 



dX 9g(0 



(CIO) 



Comparing Eqs. (C8) and (CIO), we notice that the notation p''V-yfa0 could be ambiguous. Indeed if we use the 
canonical coordinate system for the tangent space, then d/dx'^ is to be taken at p^^ fixed, but if we take the tetrad 
basis (or another coordinate system for the tangent space), then d/dx^ is to be taken with g*^*^ fixed. 

Eq. (CIO) is not exactly the desired form of the Boltzmann equation when the distribution function depends on 
(x'', q'^'^'). In fact, the use of the tetrad basis makes it natural to work with spherical coordinates in the tangent space. 
In order to introduce them, we first relate the Cartesian derivative in the tangent space, that is, the derivative with 
respect to q'^^\ to the covariant derivative on the unit sphere which is described by the possible directions n'^^^ of the 
momentum. We must stress that at any point of space-time, a;'', a distribution function (tensor- valued like /^^ or 
scalar valued like its trace /) which depends on (x^ , q^^^ ) can be considered as a field in the tangent space because 
the tangent space at a given point can be considered as a fiat three-dimensional manifold whose points are labelled 
by q^'^' and the natural covariant derivative in this manifold is d/dq^K Using that q^^' — qn^^\ and the property 






-5('') 



0)' 



it is possible to show the following relations 



8T 8T 



dS 

9oW 



dhiq 
dS 

7 

dh\q 



,W 



D(i)S , 



with i?(,)T^, = g^(,)0)5^"^/^^ 



(i) 



with D(^,^S = qS(^i) 



U) 



dS 
dqU) 



(Cll) 

(C12) 
(C13) 



where T^i,(g(*)) is a projected tensor field in the tangent space such that Tp_v{q^^^)n'' — Tp^{q^-^^)n'' = 0, and S'((7*-*^) 
is a scalar field in the tangent space. Here -D(i) is the covariant derivative on the two-sphere associated with the unit 
direction vector n^*', and it appears naturally as an induced derivative on the sphere, given that this is the surface 
orthogonal to n'^'^ (see Ref. [51? for more details on induced derivatives). We can then deduce the useful property 



qjn Oi, 



9g(*) dlnq 



Hi) +-D(i)/^,y . 



(C14) 



Given that 



dq^ 
~dX 



= q 



ding 



,(■0 



dn(*) 
dA 



(C15) 



from Eqs. (C12) and Eq. (CIO), we then find that the Boltzmann equation takes the form 

dx'^ dfp^ ding 



c aa '^Jafi _ c a q Pyi r 

^M ^^ ^^ -*P *- ^^^"^dA ■ 51ng dA 



dn(') 



-^(O^M'^^lj^ - ^M!' ■ 



(C16) 
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3. Decomposition of the Boltzmann equation 



In order to obtain equations for the components of /^i,, I^P^v and V , we want to apply the same type of decom- 
position on the equation itseh. Applying V/VX on the decomposition (2.18) of /^^ leads to 



Vf^, I (VI 
VX 2[VX^'"' 



vs^A , pp., , i 



'''-^'^'S.... + I^]+^ + -e^,.,['—e,ofn^ + V 



vv 



VX 



(C17) 



We then need to screen-project this equation in order to obtain the tensor-valued Boltzmann equation. Then the last 
two terms vanish. Indeed, first 






dA^ J- ^^^rj^pV-^vT^^vr''^ 



= 0. 



(C18) 



Second, from the normalization condition of em)^ and rt'', we can show that 



g(0) T^H^f 
^ ^ VX ~ ' 



Vn^' „ 

n,, — . 

^ VX 



(C19) 



This means that the derivative of e(o)'^ is orthogonal to e'"-*^ and the derivative of n^ is orthogonal to n^. We thus 
find that 






VX 



n + 6(0) 



,Vrf 
VX 



(C20) 



Finally, we obtain that the Boltzmann equation for the tensor-valued distribution functions can be split into the 
desired form as 



J/x "" PA ~ 2VX " VX 2^VX "^ ''^^ — L-A"' ■ 



(C21) 



4. Expression of the Boltzmann equation for polarisation 



First of all, the Boltzmann equation for the circular polarisation is the same as that for the intensity, but with 
V = Q and with a vanishing collision term since it is not generated by the Compton scattering. We shall not study 
further the equation dictating the evolution of V since it should remain null at all time unless generated by other 
types of coUisions. 

As for the linear polarisation, let us write down the basic equations for the tetrad components as commonly done 
at first order. By moving to the tetrad components, one can rewrite the covariant derivative in the Liouville operator 
as 



= A", i^rrpi c Wc [d) ( p Je) dx^ , 9P(c)(d)dlng 

3(a)^e(6) £[PJp, = i(a)' 'S(^b)' [P{c){d)\{e)e' \-^ + -g^ ^ 



d?l(^) 

P'i^)PicKd)-^ 



(C22) 



where 



/(a)(6)|(c) = e(cfdf,fi^a)ib) " w'^ ' (a)(c) f (d)(b) " W^ ' (b)(c)f(a)(d) 



M), 



(C23) 



and u'(°)(f,)(c) is the Ricci rotation coefficient defined by w^°'\i,)(^c) = g^°'\^ {c)G{b)^ ■ Here we also notice that only the 
spatial component has non-vanishing term because the projection of S'^'' onto e^'m) vanishes by construction. After 
the decomposition of the Boltzmann equation for polarisation, one obtains the following equation for the temperature 
part as in Eq. (3.25) up to the second order 



C 



ii)U) 



^r, 



WO) 



\r)(]),kn 



(fc) 



dyy 
dA 



(1) 



^W(j)' 



~dX 



(1) 



dn^*^) 



W(i),fc 



dA 



(1) 



D(k)P(i)(j) 



. (1) 



a{w^''\iW) + ^"^'^OW^^^O (%)^'^^(^)0-) + ^O)^'^^^^)^) 



(C24) 
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As for the spectral distortion, the equation is given by 

The complete expression of the collision term for the linear polarisation is given by 



(C25) 






neCJT ( -^(,)0) - 47^^)0-)^'^^'^ [{Qrn^kW)) - 2(nfe)(0>J + ^'*'^"(/c)7'(.)0) 



+ aneaT 7^.)0)^"^('' 



^(fc)(0 ^^(fe)(0 



-3C®-P(i)(j)+(5eC(^)(j) 



^WO) = «"eaT -%a) - t'7^,)(,)('=)('^ {ym^km) - 2(y(,)(,)) + 7^.)0)('^)(''Qf,.)(,) - C«7'(,o(,) 



where 7(2)0) ''^■'^''' ^^ ^ traceless projection operator with respect to S*, 



WO) 



'(«)(j) - '^(«) -^O) 2 '^(*)0) 



(C26) 

(C27) 

(C28) 



Finally before closing this section let us explicitly write down the equation for the temperature part for the sake of 
completeness. It is of the form 






and using that at first order 



the explicit form of Cf^-., •, is 



^^''\o(o) = -/?''',/] , w^^\i){r,^) = -hj'', 



(C29) 



(C30) 



- (a.fc - A,fcnW + /ife/n« + 2;i„[fc,,p«n("))i?('=)7'(,)o) 



(C31) 



Appendix D: Technical details about linear polarisation 
1. Gauge transformation of a tensor-valued distribution function 



In this section, we investigate the transformation property of the distribution matrix up to the second order. The 
transformation properties of / and V, which are scalar distribution functions, have already been investigated in the 
main text. In this section we focus on the case of the linear polarisation and set / = V^ = such that /^^ = P^^. 
Since the polarisation matrix is at least first order, this means that we need only to keep terms which are first order 
in the coordinates transformation (T, L*). 

First we must understand the transformation properties of the screen projector. As for the tetrads, it is not a 
purely geometric quantity and it is not invariant under a change of coordinates. Indeed, it is defined with respect to 
the time-like tetrad which depends on the choice of coordinates. We thus have the two related screen projectors 



S'(p)=g + e(°)®e('')-n(»n, S{p) = g + e^°^ ®e^"^ -h^fi , with n=^^^e^°\ h = ^^ - e^°K (Bl) 

pf^ey 



r^ 



If we use the transformation rule at the same point (4.12) wc find that at first order in the transformation the relation 
between these projectors is given by [55] 



V ° 5(*,'?^°^) = V(a;,<Z<''^) + 2 (e["^ +n(,,) 5,)„(a;, g('^))y" , with V = -A'^U,,, = -d'Te, 



iW ^ 

(O)^(i) 



Note that if we use tetrad coordinates for the argument of the screen projectors, then 



S(g,nW)=g + e 



(0) fs, ^(0) 



n(,)n(^)e 



(*) 



,U) 



So c{q,nS''') = g + e 



p(") 6?. ^(0) 



"W^O^e 



£(i) « ^U) 



(D2) 
(D3) 
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This means that when the coordinates of the projectors are expressed in the tetrad basis associated with the corre- 
sponding coordinates system, we obtain 



S({)(j){q,n'^''^) ^6ij -n(,)7 



■0) 



= 5; 



Wb) 



(D4) 



Thus in any gauge, the expression of the related screen projector in the tetrad basis is the same by construction, and 
S'jijQ) depends actuaUy only on n'-''\ However it must remain clear that these two tensors are geometrically different 
since they are associated with different coordinates systems and indeed their relation is given at first order by (D2). 



The distribution tensor is also not geometrically invariant since for every observer used to define the screen projector, 
we must consider a different distribution tensor. However all the possible distribution matrices are related through 
projections and we find that the distribution tensors defined by the tetrad e^p) and C/q) are related at the same point 
of the tangent bundle by [28] 

/;r o £{i, gW) = 5^ o £{x, q^^)St o c{x, q^^)flf {x, qW) . (D5) 

This means that the only requirement is to project the distribution function so that it is projected with respect to 
the new observer and the new direction. Combining this transformation rule with (D2) we obtain at first order the 
transformation rule as 



M 



/;<"' o c{x, gW) ^ ;;<") (^, ^W) + 2 [e^;i + n^A /;<:> [x, gW)^/ 



=(") 



M\\ra 



(D6) 



6(0) 



If we project this expression onto the tetrad components, and noting f(i)(j) = //-/u -i and f{i)(j) = fa)/.-): we obtain 



'WO)' 



find the transformation rule under a gauge transformation, we need to expand the left hand side around {x, q^^'). At 
first order we get 



/^(~^o5(*,gW) 



i + ^^a^ + ^'^'^a|^)%)(7)°^(^''^'^) 



/(^)(j^o£(a;,gW) 



^''a^+^^'^a|^)^«(^) 



(D8) 



Using the expansion (C12) we then obtain the gauge transformation rule for the tensor valued distribution function 
in tetrad coordinates. Noting f(i)(j) = f7~)p) ° c for simplicity, this reads 



d 



/w(j) = /wo) - U''/):;:^: + '^i'^«7)i^ + '^"^*^^w ) •^wo") ^ /ww-^ m^u) ^ hm)^ ,'ij^w > 



dh\q 



-[k Ml) 



-[k Ml) 



(D9) 



where it is implied that all quantities are evaluated either at x or at {x,q,n^'''>). 

For completeness we report the explicit form of the gauge transformation for 'P(i)(j) which is obtained from the 
above transformation rule and the spectral decomposition (2.261 



V^ 



h)U) = ■^(^)O) - ^^'',1] {s^^\i)'^(k)U) + S^^\j)'J^{k){t)) - Cdf,V(^i)^j) - (5n('=)L>(fc)7'(^)(j) 



(DIO) 



It is also found, as expected, that the spectral distortion ^(mj) part is gauge invariant since it vanishes on the 
background and at first order. 



2. Gauge transformation for Liouville and Collision terms 



We deduce from the transformation rule (DIO) and the spectral decomposition (3.22 1 and (3.24) that the spectral 
distortion part, 'C^^wx and C,^w -., must be gauge invariant. Concerning the temperature part, they should transform 



as (noting C^^^^,^^ = C^^^^ o 5 and Cf^^^^^ 



^ii)U) ^ ^(i)U) ^ ^'''■'l ('^^''W'^W(J) + ^'^^\j)^{k){i)) - ^^^M^W(j) - '5"*''-'^(fe)'C(J)(j) 



^mj) = ^mj) - ^'^'1 {^^^\^)^fk)ij) + s^^\j)Ci'k)(^)) - ^""d^^mj) - '5"^''^^(fc)Cw(j) 



(Dll) 
(D12) 



Using the transformation rules derived in this paper, we checked that this is indeed the case when using the detailed 
form of the Liouville and collision operators. 
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Appendix E: Extraction of temperature and spectral distortion 

The functions y and Q can be extracted thanks to the integrals of the type 

//q2+»dg 



Mn[f] ^ 



(3 + n) / I{q)q^+"dq ' 



just by applying them order by order to I{q), using that AiolV^I] = 0. We then obtain 



i2;(2)=l(M,[/(2)]-M„[/(2)])-leW = 



(El) 

(E2a) 
(E2b) 

(E2c) 



Similarly to what can be done for the intensity part, the spectral components of polarisation can be extracted 
thanks to 






(E3a) 
(E3b) 

(E3c) 
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